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ABSTRACT 



A boundary layer analysis has been presented to study the combined effects of magnetic field, viscous dissipation 
and heat source/sink on the fluid flow and heat transfer of a power law fluid flowing over a wedge, with magnetic field is 
investigated. The governing partial differential equations of the hydro magnetic free convective boundary layer flow are 
reduced to ordinary differential equations by the applications of group theory. By using Quasi -linearization technique first 
linearized the coupled non-linear equations, and then they are solved by numerically by an implicit finite difference 
method. Numerical solutions for the governing momentum and energy equations are obtained. Results are presented in 
terms of the velocity profiles and temperature profiles for the different flow parameters, such as, Magnetic field parameter 
M, Prandtl number Pr, Eckert number Ec, the flow behavior index n, and the wedge angle parameter m and heat source 
parameter S. Variation of heat transfer and skin friction for different values of Ec, Pr, M and m are presented. Heat transfer 
and skin friction results are compared for various values of flow behavior index n governing the nature of the fluid and also 
for different wedge angles. 

KEYWORDS: Non-Newtonian, Flow Behavior Index, MHD, Eckert Number, Heat Source, Finite Difference Scheme 
INTRODUCTION 

Non-Newtonian fluid is a fluid in which the viscosity changes with applied strain rate. As a result, non -Newtonian 
fluids may not have well defined viscosity. In modern technology and in industrial application, non -Newtonian fluids play 
an important role. Many processes in modern technology use non-Newtonian fluids as working fluids in heat exchanger. 

The prediction of heat transfer characteristics of these fluids have been studied widely during the past decades. 
Due to the growing use of these non-Newtonian substances in various manufacturing and processing industries. During the 
past 50-60 years, there has been a growing recognition of the fact that many substances of industrial significance, 
especially of multi-phase nature (such as foams, emulsions, dispersions and suspensions, slurries) and polymeric melts and 
solutions (both natural and manmade) do not conform to the Newtonian postulate of the linear relationship between shear 
stress and rate of shear in simple shear. Accordingly these fluids are variously known as non-Newtonian, on-linear, 
complex or rheologically complex fluids. Indeed so widespread is the non -Newtonian fluid behavior in nature and in 
technology that it would be no exaggeration to say that the Newtonian fluid behavior is an exception rather that the rule. 
Since 1960, a considerable attention has been devoted to predict the theoretical analysis of an external boundary layer flow 
for power law non-Newtonian fluids were first performed by Schowalter[l] and Acrivos et al [2] .The study of heat and 
mass transfer in a non - Newtonian power law fluid obeying the Ostwald-de Waele rheological model, 



16 



P. Kavitha & N. Kishan 




3v. 



By 



has been attracting the interest of researchers and scientist in the recent time due to its applications in food, 
polymer, petro-chemical, geothermal, rubber, paint and biological industries. The two parameter rheological equation (1) is 
also known as the power law model. When n=l, the equation represents a Newtonian fluid with a dynamic coefficient of 
viscosity m. Therefore, deviation of n from a unity indicates the degree of deviation from Newtonian behavior. For n < 1, 
the fluid is pseudo plastic and for n > 1, the fluid is dilatant .n is power law exponent and m is the consistency coefficient. 
The flow of a fluid past a wedge is of fundamental importance since this type of flow constitutes a general and wide class 
of flows in which the free stream velocity is proportional to a power of the length of coordinate measured from the 
stagnation point. In these types of problems, the well known Falkner-Skan transformation is used to reduce boundary- 
layer equations into ordinary differential equations for similar flows. It can also be used for non -similar flows for 
convenience in numerical work because it reduces, even if it does not eliminate dependence on the x-coordinates. The 
solutions of the Falkner-Skan equations are sometimes referred to as wedge-flow solutions with only two of the wedge 
flows being common in practice. The dimensionless wedge angle parameter plays an important role in such type of 
problems because it denotes the shape factor of the velocity profiles. Acrivos et al [2] considered the problem of natural 
convection heat transfer to power-law fluids for different geometric configuration. The application of boundary layer flow 
to power law pseudo-plastic fluids, similar solutions has been studied by Scho waiter [l].The flow of a power law fluid past 
a symmetrical wedge is studied by S.S.Elghabaty and Abdel-Rehman [3] in the neighborhood of the stagnation point when 
an external magnetic field is applied. Kuo Bor-Lih [4] studied the effects of heat transfer analysis for the Falkner-Skan 
wedge flow by the differential transformation method. W.T.Cheng and H.T.Lin [5] analyzed the non- similarity solution 
and correlation of transient heat transfer in laminar boundary layer flow over a wedge. The transient thermal response of a 
power-law type non-Newtonian, laminar boundary flow over a wedge is investigated by Rama Subba Reddy [6]. The flow 
and heat transfer characteristics of a second-order fluid over a vertical wedge with buoyancy forces have been analyzed by 
M.Kumari, et al [7]. For the non-Newtonian power-law fluids, the hydro magnetic problem of the MHD boundary layer 
flow over a continuously moving surface over a stretching sheet has been dealt with by several authors e.g. Anderson et al 
[8], Cortell [9] and Mahmoud [10]. The effect of magnetic field is found to decrease the velocity distribution and thus to 
increase the skin-friction coefficient. However, relatively less attention has been paid to the accompanying heat transfer 
problem of power-law fluids past a stretching surface in the presence of a magnetic field. The MHD free-convection flow 
of non-Newtonian power-law fluid at a stretching surface with a uniform fee-stream was studied by Emad M. Abp-Eldahab 
and Ahmed M. Salem [11]. 

The non-Darcy flow characteristics of power-law non-Newtonian fluids past a wedge embedded in a porous 
medium have been studied by Youn J.Kim [12]. Ramamurthy [13] studied convective heat transfer to power law, 
non-Newtonian fluids flowing over a wedge with viscous dissipation. MHD laminar boundary layer flow over a wedge 
with suction or injection had been discussed by N.G.Kafoussias and N.D.Nanousis [14] and M.Kumari [15] discussed the 
effect of large blowing rates on the steady laminar incompressible electrically conducting fluid over an infinite wedge with 
a magnetic field applied to the wedge. A non similar boundary layer is presented of mixed convection in power-law type 
non-Newtonian fluids was studied by Rama subba reddy [16], a boundary layer analysis has been presented by 
M.A.Mansour et.al[17] Mehmet Akcay[18] studied Flow of power-law fluids over a moving wedge surface. 
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In the present work, we investigated the effect of MHD heat transfer to non-Newtonian power-law fluids flowing 
over a wedge with viscous dissipation. The flow is subjected to external magnetic field. The governing flow equations have 
been solved by Quasi-linearization technique with finite difference scheme. 

MATHEMATICAL ANALYSIS 

A steady, two dimensional incompressible laminar boundary layer flow of a non -Newtonian fluid, which obeys 
the power law model, past a wedge of included angle n P is shown in figure 1. The wedge is maintained at a uniform 

surface temperature t w and the free stream fluid has a temperature t x . The governing equations of continuity, momentum 

and energy are: 




Figure 1 
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Write, the dimensionless temperature q _ t — t^ 

t,., -t„ 



And the shear stress T = n (— )" 

dy 



(1) 
(2) 
(3) 

(4) 
(5) 



In (5),ni is called the consistency index and n the flow behavior index. For shear thickening dilatants fluids n has 
a value greater than 1 (n>l) and for pseudo-plastic fluids n value less than l(n<l). 



The boundary conditions are 
At y = 0, u = v = and 0=1 
As,y— >oo u = and 0=0 
The free stream velocity is taken as = Cx m 
Where, the wedge angle parameter _ r / 

Transformation using group theory 

Let x= al x,y = a 1 p y -Independent variables u = a[u, v = a[v,U = a™U ^ - Dependent variables 
On substitution of the above, the continuity and moment equations will be 



(6a) 
(6b) 
(7a) 

(7b) 
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^^^^^^^H^i,^^.^^; (9) 

dx dy p dy dy 2 dx p 

To establish variance it is necessary to choose the exponents of the parameter in such a way that the 
transformed equations look exactly the same in terms of the new variables. For achieving this, it is necessary to equate the 
exponent and cancel out the common powers of ai.This gives 

1 - r = p - s From which 

1 — 2r = p — s — r = p ( 1 + ri) — nr = 1 — 2 vv 

pO +//) / (10) 

/■ = vv ' = — 

Oz-2) 



/7( 2// 1 ) /( // 1 ) 



(ID 



It is now necessary to find an absolute invariant rj 1 which does not vary under the transformation of the 
independent variable. 

Let n x = yx h = yx h a~ p ~ lb =yx^ (v-p-lb = 0) (12) 

The new dependent variables to replace u and v are now sought and they are to be such that they are invariant 
under the complete set of transformations. Choosing 

— r 

d -r-ld d ~T 

g x = ux = a x u X = ux 1 

— s 

e —s—le — e r 

g 2 = VX = a x V X = VX 1 

g 3 U^x h =a^ w ~ lh U x h = U^x~ 

And considering gi9 an d g 3 being functions of rj x , the partial differential equations can be reduced to a set 
of ordinary differential equations. 

Let gi = a^X g 2 = Hir, x \ g 3 =/(//!) then 

u = Gir/^x 1 (13) 
v = H^r/^x 1 (14) 

And Uao = nn^x* (15) 
from (7a) and (15) 

/(t 7i ) = C and = fji (16) 
/ 

Substituting (13), (14) along with (12) in continuity equation (1) we get 

XL I 

Primes denote differentiation with respect to rj l 
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From (10) and (11) 



= x 

x 



s-p 

I __ v Hn-T) 



Hence, the continuity equation will be _ p r (18) 
After substituting the transformed variables in the momentum equation (2) we get 

x T ~ 1} [G(y - ^ r, x G> )] + [HG< ] = 

r-p r-2p ,2w D 2 r (19) 

^[^— ]-^— [G']-G" + C 2 x ( -- 1) ^-^G^ 
P / p 

On substitution of the values for r and s from (10) and (11) the exponents of x terms in equation (19) are found to 
be equal, that is 

2r_ x s+r-p rn-p(\+n) 2w^ 2p(l+n)-ln 

x T ~ = = x ' = x T ~ =x ' ( '- 2) (20) 

The momentum equation will then be 

11 p I p 

On perusal of equation (20), it can be seen that 

r xv p(\ + n) I (22) 

11 l(n- 2) 

p_ m(n - 2) + 1 (23) 
/ (1+//) 

And s_ = "'(2/7 -Y)-n (24) 
/ (1 + n) 

Equations (22) and (23) in (18) will change the continuity equation to the form: 

^ t= m(n-2) + l G (25) 

(1 + ri) 

Let G = F\r ll \ G'^F'^), G" = F"^) (26) 

The equation (25) can be written as 

m O - 2) + 1 

// = / ;, F —i nF 

(1 + n) 

On integration 

/7 = m(n-2) + l F , m(2/7-l) + l (27) 

(1 + 77) (1 + /7) 



Substituting (22), (23), (26) and (27)in (21), the momentum equation will be 
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m(^) 2 - m ^r~^ + 1 ^"= (28) 
(1 + rc) 

Hence, it is seen that the continuity and momentum equations combine and reduce to equation (28). Defining new 
variables/; = d 1 rj 1 where 

n y < 29 > 

'1 m(n-2)+l 

^(Vi) = KfW>F'<nJ = KdJ'irj) (30a) 

= k^fitj) ( F"'(^) = Mi7'"07) (30b) 
On substituting (30) in (28), the momentum equation will be as follows : 

mkX(f) 2 ~ m(2 " "*? + 1 *,W= ^(M. 2 )" 4 M. 3 (/"r/ , "+mC 2 - ^ W*^ < 31 > 
(1 + n) p p 

Choosing = C (32) 



<2, =[^C 2 -"] ( 



(33) 



And similarity variable 

„ Z ^ = ^ (34) 



( 



Where ^ = 2?^ jc +/ Equation (31) reduces to 

^_ (/ , y + m (2n-l) + l # „ +m[1 _ (/ , )2] _ M/ , = (35) 

The primes denote differentiation with respect to rj. From (13), (22), (26), (30) and (7a) the velocity component u 
can be written as: 

u = cx m r(7 7 ) = u O0 r(? 7 ) (36) 

From (14), (24) and (27), the velocity component v can be written as 

v = C [ m(/i-2) + l m(2^-l) + l /]x OT( ^ (37) 
d l (1 + rc) (1 + rc) 

Using the values of u and v from (36) and (37), the energy equation (3) can be reduced as 

"+A 1 .Pr./79' + (f"y +n .Ec.Pr-hSPr0 = O (38) 
Where primes denote the differentiation with respect to r| and 
Modified Prandtl number _ r n , u^x 



Pr = [ 1 ] 
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Eckert number _ u _ 2 

Ec 



Heatsource/sink s _ Q(t^ — Q and 

A _ /7/(2a7 — 1) + 1 

The transformed boundary conditions are 

At ^ = o,y = /' = 0,<9 = 1 (39a) 

As ry > oo,/' = 1,<9 = O (39b) 

Quasi-linearization of all highly non-linear terms Bellman & Kalaba [19] in the equation (35) gives 

APf ' ' '+Aiy ' -(2mF f +M )/'+A 1 F' ' / = A x FF"-m(\ + (F' ) 2 ) (40) 

Where = n(f"Y~ l A _ ^(2^ - 1) + 1 and a/? = (/ Here F is assumed to be known function, i.e., F is 

' 1_ (1 + n) 

the n th approximation to the solution and / is the (n+l) th approximation to the solution. Using Crank-Nicolson implicit 
finite difference scheme, the equations (40) and (38) are transformed to 

avvu + 1] + bmm + cuuu - u - auvu -2] = Em 

a^iWU + 1] + b^iWU] + c^iWU - 1] = -^M ( 41 ) 

= AP[/] + h\F\i\ - 0.5/z 2 (2mF'[/] + Af ), 
b\i\ = -3AP[/] + ^A^T*'] - 2AJiF[i\ 
c[i] = 3AP[/] + A x hF\i\ + 0.5/z 2 (2mF'[/] + Af ), 
£/[/] = AP[/],^H = ^(A^C/]^"^'] -m-m((^'[i]) 2 ) 

^[i] = 1 + 0.5/zA, Pr/^m, ^[z] (-2) + /z 2 Pr,S (42) 

Here h represents the mesh size taken as h =0.005. The system of equations (41) and (42) has been solved by 
Gauss-Seidel iteration method and numerical values are carried out after executing the computer program for it. In order to 
prove convergence of finite-difference scheme, the computations are carried out for different values of h by running same 
programme. Negligible change is observed, after each cycle of iteration the tolerance set asl0" 4 is satisfied at all points. 

Heat Transfer and Skin Friction 

Heat transfer and skin friction values are very important from engineering point of view. The local heat flux at the 
surface can be calculated from the Fourier's Law. 

q w =-K^\ y=0=Kt w -U (43) 
dy 

Here primes denote the rate of heat transfer per unit area. The local Nusselt number if given by 

CAfE/) = IhL (44) 

X k 

In terms of transformed variable 
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(NU) X = [ P^ 2 "^ mC2 " W)+W ] i + n [_*•] (45) 
or (Aft/), _ (46) 

— - J 77=0 

(Re) x \+n 

The shear stress at the wall is given by 

r ,a^, | I* (47) 

And in terms of transformed variables 

This can be expressed as skin friction coefficient by the following equation 



_ 2t w or c / (Re) x i+« - r/-n« 

C r — ' 

f pUl 2.0 

RESULTS AND DISCUSSIONS 

The system of non-linear equations (35) and (38) are solved numerically under the conditions given by (39), using 
a marching technique and by applying the Crank-Nicolson implicit finite difference method. The resulting system of 
difference equations has to be solved in the infinite domain < r\ < oc. A finite domain in the iq- direction can be used 
instead with v\ chosen large enough to ensure that the solutions are not affected by imposing the asymptotic conditions at a 
finite distance. The independence of the results from the length of the finite domain and the grid density was ensured and 
successfully checked by various trial and error numerical experiments. Computations are carried out for r^lO which is 
found adequate for the ranges of the parameters studied here. In order to validate our method, we have compared our 
results for the Newtonian case (n=l) with those of Ramamurthy et.al [13] and found excellent agreement. Details of the 
velocities, temperature fields were presented in curves for various values of the different parameters of the problem, 
e.g., the flow behavior index n, magnetic field parameter M, wedge angel parameter m, Prandtl number Pr, heat source 
parameter S and Eckert number Ec. 

Figure 1(a) - (b) respectively, depict the effect of shear thinning (n<l), Newtonian (n=l), and shear thickening 
(n>l) fluids on the horizontal velocity profiles /' with r|. From the geometrical representation, we notice that increasing 

the values of the flow behavior index parameter n is to reduce the velocity profile / ' , which tends to 1 as the space 
variable r| increases from the boundary surface, in the case of (n<l). Further, it is observed that the velocity profiles 
f r increases with the increase of n in the case of n>l . 

Figure 2(a)-(b) presents the velocity profile /' for (a) m=1.5 and (b) m=0.6 for various values of magnetic 
parameter M. These figures indicate that an increase in M results in an increase in the velocity profiles /' . However the 

effect of M is more when m=0.6 in comparison to its effects when m=1.5. The graphs for the velocity profiles /'with r| for 
n=l for different values of wedge angle parameter m are shown graphically in figure-3. The effect of increasing value of 
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wedge angle parameter m is to increase the velocity profile. Figure. 4 illustrate the dimensionless temperature profiles for 
different values of Prandtl number for Ec = 0.1 and (a) n = 0.8 (b) n=l(c) n=1.6 in the presence of constant wedge angle 
parameter. The temperature of the fluid decreases with increase of Prandtl number Pr for different values of flow behavior 
index n. Figures. 5(a)-(c) represent the dimensionless temperature profiles for different values of Prandtl number, for 
Ec=1.0 and various values of wedge angle parameter m. The temperature of the fluid decreases with increase of Prandtl 
number Pr is observed. Figure 6 (a)-(c) shows the dimensionless temperature profiles for different values of Prandtl 
number for Ec = 10 and it is seen that the temperature of the fluid increases with the increase in Prandtl number. Figure 7 
(a)-(c) portrays the dimensionless temperature profiles for different values of flow behavior index n, when m = 1.0 and Pr 
= 1.0. for different values of Ec=0.1, Ec=l, Ec=10 is shown respectively. It is noticed that the temperature of the fluid 
decreases with increase of flow behavior index n. This effect is more when Ec=10. Figure 8 (a)-(c) represent the 
dimensionless temperature profiles for different values of Eckert number for n=0.8, n=l, n=1.6 respectively. The 
temperature of the fluid increases with increase of Eckert number Ec noticed. From figure 9 decrease in temperature of the 
fluid is noticed with the increase of wedge angle parameter m for fixed values of n and Ec. 

Figure 10(a)-(c) represent the temperature profiles for different values of heat source/sink S for n=0.8, n=l and 
n=1.6 respectively. The effect of heat source/sink parameter S on temperature is executed. From figure it is noticed that the 
dimensionless temperature decreases for increasing strength of heat source and due to increase of heat sink strength the 
temperature increases. Figure 11 reflects the skin friction profiles C/ ( Re )% + i for different values of wedge angle 

parameter m. It is observed that the skin friction increases with increase in wedge angle parameter. Figure 12 represents the 
skin friction profiles c (Re)^" +1 for different values of magnetic field parameter M. It is noticed that the skin friction 

increases with increase of M. The variation of heat transfer coefficient (mjj} f° r different values of wedge angle 

(Re) X + i 

parameter m is shown in Figure (13). It is observed that the skin friction increases with increase of wedge angle parameter 
m. The heat transfer is enhanced for the increase of magnetic field parameter M is observed in Figure (14). 




Figure 1(a): Velocity Profiles for Figure 1(b): Velocity Profiles for Figure 2(a): Velocity Profiles for 
Different Values of Different Values of Different Values of M 

m=l, Pr=l, Ec=l, M=0, S=0 n, m=l, Pr=l, Ec=l, M=0, S=0 n=l, Pr=l, Ec=l, m=1.5, S=0 
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Figure 2(b): Velocity Profiles for 
Different Values ofM 
n=l, Pr=l, Ec=l, m=0.6, S=0.0 



Figure 3(a): Velocity Profiles for Figure 4(a): Temperature Profiles for 
Different Values of m for Different Values of Pr 

11=1, Pr=l, Ec=l, M=0, S=0.0 n=0.8, m=0.1, Ec=0.1, M=0, S=0.0 
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Figure 4(b): Temperature Profiles for Figure 4(c): Temperature Profiles for Figure 5(a): Temperature Profiles for 
Different Values of Pr for Different Values of Pr for Different Values of Pr for 

n=l, m=0.1, Ec=0.1, M=0, S=0.0 n=1.6, m=0.1, Ec=0.1, M=0, S=0.0 n=0.8, m=0.1, Ec=l, M=0, S=0.0 
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Figure 5(b): Temperature Profiles for Figure 5(c): Temperature Profiles for Figure 6(a): Temperature Profiles for 
Different Values of Pr for Different Values of Pr for Different Values of Pr for 

n=l, m=0.1, Ec=l, M=0, S=0.0 n=1.6, m=0.1, Ec=l, M=0, S=0 n=0.8, m=0.1, Ec=10, M=0, S=0 
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Figure 6(b): Temperature Profiles for Figure 6(c): Temperature profiles for Figure 7(a): Temperature Profiles for 
Different Values of Pr Different Values of Pr Different Values of n 

n=l, m=0.1, Ec=10, M=0.0, S=0.0 n=1.6, m=0.1, Ec=10, M=0, S=0.0 m=l, Pr=l, Ec=0.1, M=0, S=0.0 



n=0.8 

n=l 

n=1.6 





2 
o 

1 



Ec=0.1 
• Ec=l 
Ec=10 



2 4 6 

71 



Figure 7(b): Temperature Profiles for Figure 7(c): Temperature Profiles for Figure 8(a): Temperature Profiles for 
Different Values of n for Different Values of n for Different Values of Ec for 

m=l, Pr=l, Ec=l, M=0, S=0.0 m=l, Pr=l, Ec=10, M=0, S=0.0 n-0.8, m=l, Pr=l, M=0, S=0.0 
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Figure 8(b): Temperature Profiles for Figure 8(c): Temperature Profiles for Figure 9: Temperature Profiles for 
Different Values of Ec Different Values of Ec for Different Values of m for 

for n=l, m=l, Pr=l, M=0, S=0.0 n-1.6, m=l, Pr=l, M=0, S=0.0 Ec=0.1, Pr=l, n=l, M=0, S=0.0 
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Figure 10(a): Temperature Profiles for Figure 10(b): Temperature Profiles for Figure 10(c): Temperature Profiles for 
* A c Different Values of S ,A * A ^ 
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Figure 11: Skin Friction Values for Different Values of 
m for n=l, Ec=0.1, Pr=l, M=0, S=0.0 



Figure 12: Skin Friction Values for Different Values of M for 
n=l, Ec=0.1, Pr=l, m=0.1, S=0.0 
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Figure 13: Heat Transfer Coefficient for Different Figure 14: Heat Transfer Coefficient for Different Values of M 
Values of mforn=l, Pr=l, Ec=0.1, M=0.0 for Pr=l, Ec=0.1, m=0.0 

CONCLUSIONS 

The problem of magneto -hydrodynamic flow and heat transfer to non-Newtonian power-law fluid past with heat 
source/sink in the presence of viscous dissipation is studied theoretically. The governing partial differential equations are 
transformed into ordinary differential equations by using appropriate similarity transformation and resulting boundary 
value problem is solved numerically by second order finite difference scheme. The effects of various governing parameters 
such as the Magnetic field parameter M, Prandtl number Pr, Eckert number Ec, the flow behavior index n, and the wedge 
angle parameter m and heat source/sink parameter S. Variation of heat transfer and skin friction for different values of Ec, 
Pr, M and m were examined. The following conclusions are drawn from the computed numerical values. 



Quasilinearization Approach to MHD Heat Transfer to Non-Newtonian Power-Law Fluids Flowing Over a 27 
Wedge with Heat Source/Sink in the Presence of Viscous Dissipation 

• The heat transfer coefficient value reduces as the power-law index parameter increases. 

• The skin friction value decreases as the power-law index parameter value increases. 

• The effect of magnetic field decreases the skin friction values. 

• In dilatant fluids the skin friction and heat transfer rapidly decrease as n increases. 

• Temperature decreases for increasing strength of heat source and due to increase of heat sink strength the 
temperature increases. 
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